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ABSTRACT 

We review the properties of young superstellar clusters and the impact that 
their evolution has in their host galaxies. In particular we look at the two different 
(Nj [ star-forming feedback modes: positive and negative feedback. The development 



of strong isotropic winds emanating from massive clusters, capable of disrupting 



the remains of the parental cloud as well as causing the large-scale restructuring of 
the surrounding ISM, has usually been taken as a negative feedback agent. Here 
we show the impact that radiative cooling has on the resultant outflows and then, 
as an extreme example, we infer from the observations of M82 the detailed inner 
structure of supergalactic winds and define through numerical simulations the 
ingredients required to match such structures. We also show how when radiative 
cooling becomes significant within the star cluster volume itself (~ 30% of the 
deposited energy), the force of gravity takes over and drives in situ all the matter 
deposited by winds and supernovae into several generations of star formation. A 
situation in which the mass deposition rate Msc, instead of causing a wind as in 
the adiabatic solution, turns into a positive feedback star-forming mode equal to 
the star formation rate. 



1. The properties of superstellar clusters 

The discovery by HST of a large population of unusually compact young superstellar 
clusters (SSCs) within starburst galaxies (see review by Ho 1997; Johnson et al. 2001; Colina 
et al. 2002; Larsen & Richtler 2000, Kobulnicky & Johnson 1999 and the proceedings edited 
by Lamers et al. 2004), has led us to infer the unit of massive or violent star formation. 
SSCs with masses in the range of a few xlO 5 M to up to 6 xlO 7 M Q (see Walcher et al. 
2004; Pasquali et al. 2004) within a small volume of radius 3 to 10 pc, are indeed some of the 
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most energetic entities found now in a large variety of galaxies. Note also that collections 
of them have now been found within a single starburst galaxy, as in M82 and the antennae 
galaxy (see Melo et al 2005; Whitmore et al. 1999). 

The potential impact that these new units of star formation may have on the ISM of 
their host galaxies has been inferred from cluster synthesis models (see e.g. Cervino & Mas- 
Hesse 1994, Leitherer et al 1999). A coeval cluster with 10 6 M in stars, an initial mass 
function (IMF) similar to that proposed by Salpeter and an upper and lower mass range for 
the coeval event between 100 M and 1 M , leads initially to several tens of thousands of O 
stars. These however begin to disappear rather quickly (after t ~ 3 Myr) as they complete 
their evolution and explode as supernovae (SNe). The cluster evolution is so rapid that after 
10 Myr there are no O stars left within the cluster. All massive stars undergo strong stellar 
winds and all of them with a mass larger than 8 M Q will end their evolution exploding as 
supernova. Thus, one is to expect from our hypothetical cluster several tens of thousands of 
SNe over a time span of some 40 Myr. During the supernova phase a 10 6 M stellar cluster 
will produce an almost constant energy input rate of the order of 10 40 erg s -1 . On the other 
hand, the ionizing luminosity emanating from the cluster would reach a constant value of 
10 53 photons s _1 during the first 3.5 Myr of evolution to then drastically drop (as t~ 5 ) as 
the most massive members of the association explode as supernova. The rapid drop in the 
ionizing photon flux implies that after 10 Myr of evolution, the UV photon output would 
have fallen by more than two orders of magnitude from its initial value and the HII region 
that they may have originally produced would have drastically reduced its dimensions. Thus 
the HII region lifetime is restricted to the first 10 Myr of the evolution and is much shorter 
than the supernova phase. It is important to realize that only 10% of the stellar mass goes 
into stars with a mass larger than 10 M Q , however, it is this 10% the one that causes all the 
energetics from the starburst. Being massive, although smaller in number, massive stars also 
reinsert into the ISM, through their winds and SN explosions, almost 40% of the starburst 
total original mass. And thus from a starburst with an initial mass of 10 6 M one has to 
expect a total of almost 4 xlO 5 M violently injected back into the ISM, during the 4 xlO 7 
years that the SN phase lasts. From these, almost 40,000 M will be in oxygen ions and less 
than 1000 M in iron (see Silich et al. 2001). 

One of the features of the stellar synthesis models regarding the energetics of coeval star 
clusters is that they fortunately scale linearly with the mass of the starburst. It is therefore 
simple to derive the properties of starbursts of different masses, for as long as they present 
the IMF, metallicity and stars in the same mass range considered by the models. 

When dealing with the outflows generated by star clusters, another important intrinsic 
property is the metallicity of their ejected matter. This is a strongly varying function of time, 
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bound by the yields from massive stars and their evolution time. Thus, once the cluster IMF 
and the stellar mass limits are defined, the resultant metallicity of the ejected material is 
an invariant curve, independent of the cluster mass. Here we consider coeval clusters with 
a Salpeter IMF, and stars between 100 M Q and 1 M , as well as the evolutionary tracks 
with rotation of Maynet & Maeder (2002) and an instantaneous mixing of the recently 
processed metals with the stellar envelopes of the progenitors (see Silich et al. 2001 and 
Tenorio-Tagle et al. 2003, for an explicit description of the calculations). This leads to 
metallicity values (using oxygen as tracer) that rapidly reach 14 Z Q (see Figure 1), and 
although steadily decaying afterwards, the metallicity remains above solar values for a good 
deal of the evolution (for more than 20 Myr), to then fall to the original metallicity of 
the parental cloud, once the oxygen yield has been delivered. One of the main effects of an 
enhanced metallicity of the ejecta is to boost its radiative cooling and in such a case, massive 
clusters may inevitably enter into the strong cooling regime, to then find their stationary 
superwinds totally inhibited (see section 4). 



2. Feedback from superstellar clusters 

The close spacing between sources within a super-star cluster warrants a very efficient 
thermalization of all their winds and supernova explosions, leading to the high central over- 
pressure that is to drive both a superbubble and in some cases a supergalactic wind (SGW). 
The outflow from the star cluster surface is fully defined by three quantities: The mass and 
mechanical energy deposition rates (hereafter M$c and Esc) an d the radius that encom- 
passes the newly born sources (Rsc)- 

In the adiabatic case, the total mass and energy deposition rates define the temperature 
and thus the sound speed esc at the cluster surface. 



0.3fx E sc m 

where /i is the mean mass per particle and k the Boltzmann constant. On the other hand, 
the density of matter streaming out of Rsc is: 



_ M sc 

p -4nR sc c sc ' {) 

Thus at R S c (see Chevalier & Clegg 1985; hereafter CC85), the ratio of thermal and 
kinetic energy flux to the total flux is 



-4- 



20 r 




Fig. 1. — The metallicity of the matter ejected by coeval bursts of star formation. The 
metallicity (in solar units) of the matter reinserted into the ISM by SCs is plotted as a 
function of the evolution time. The curve is derived from the metal yields of Meynet & 
Maeder (2002) stellar evolution models with stellar rotation, using oxygen as a tracer. The 
estimate assumes a Salpeter IMF and 100 M Q and 1 M upper and lower mass limits. The 
model also assumes an instantaneous mixing between the newly processed metals and the 
envelopes of the progenitors. Note that under these assumptions the curve is independent 
of the cluster mass. 
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There is however a rapid evolution as matter streams away from the central starburst. After 
crossing r = R sc the gas is immediately accelerated across the steep pressure gradients and 
rapidly reaches its terminal velocity (V^, ~ 2csc)- This is due to a fast conversion of thermal 
energy, into kinetic energy of the resultant wind. 

In a recent communication (Silich et al. 2003, 2004), we have revised the properties of 
SSCs by solving the flow equations without the assumption of an adiabatic flow made by 
Chevalier & Clegg (1985). In this case, the steady-state solution results from solving 



(5) 
(6) 

q e ~ Q, (7) 

where q e and q m are the energy and mass deposition rates per unit volume (q e = q m = 
if r exceeds Rsc), Q is the cooling rate (Q = n 2 A) where n is the wind number density 
and A is the cooling function (a function of T and metallicity). Central values of density 
and temperature are found by iteration, with the restriction that the flow velocity ought 
to increase from km s _1 at the cluster center to the sound speed when it reaches the 
cluster surface. If this happens then the stationary solution is fully warranted. A solution in 
which what is put in (Msc) is & t all times equal to the outflow from the star cluster surface 
(4 tiR sc Psc c sc)- Radiative cooling leaves almost unaffected the wind density and velocity 
distributions but its temperature may for energetic clusters fall rapidly to T ~ 10 4 K, close to 
the SSC surface. Figure 2 compares the adiabatic (dotted lines) with the radiative solution 
(solid lines) for a cluster with a E sc = 3 xlO 41 erg s _1 , in which the wind temperature 
instead of steadily falling as r -4 / 3 , it rapidly falls to T ~ 10 4 K, reducing drastically the size 
of the resultant X-ray emitting volume. The energetic winds lead to a total restructuring 
of the surrounding ISM, and in some extreme cases their energetics may even reach the 
outskirts of a galaxy causing a galactic wind. 



1 d 

r 2 dr 



(p-ur 2 ) 



du 



dP 



q m u 



1_ d_ 

r 2 dr 



pur 



u 7 



P 



1 O 



- 6 - 




2- 

I , , , I , , , I , , , I , , , I , , , 

20 40 60 80 100 

R(po) 



Fig. 2. — Adiabatic vs radiative superwinds. Density, velocity and temperature distributions 
of the wind produced by a 10 pc radius SSC that deposits 3 xlO 41 erg s -1 under the adiabatic 
(dotted lines) and radiative (solid lines) assumptions. 
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We have also shown (see below section 4) that in the Esc vs star cluster radius diagram, 
there is a threshold limit that massive and energetic compact clusters may cross to find out 
that radiative cooling inhibits their stationary outflow condition (Silich et al., 2004) and 
then the matter ejected by the stellar sources, unable to escape, is to accumulate within the 
star cluster volume (see section 4). 

3. Negative feedback and the physics of supergalactic winds 

In all cases below the threshold line, which may be considered quasi-adiabatic or strongly 
radiative, the energy dumped by the central starburst, is to cause a major impact on the 
surrounding gas. The supersonic stream leads immediately to a leading shock able to heat, 
accelerate and sweep all the overtaken material into a fast expanding shell. In this way, as 
the whole structure grows, the density, temperature and thermal pressure of the wind drops 
as r~ 2 , r~ 4 / 3 and r ~ 10 / 3 , respectively (CC85). 

Note however that such a flow is exposed to the appearance of reverse shocks whenever 
it meets an obstacle cloud or when its thermal pressure becomes lower than that of the 
surrounding gas, as it is the case in strongly radiative winds and within superbubbles. There, 
the high pressure acquired by the swept up ISM becomes larger than that of the freely 
expanding ejecta (the free wind region; FWR), where p, T and P are rapidly falling. The 
situation leads to the development of a reverse shock and with it to the thermalization of 
the wind kinetic energy, reducing the size of the free-wind region. Thus for the FWR to 
extend up to large distances away from the host galaxy, the shocks would have had to evolve 
and displace all the ISM, leading to a free path into the intergalactic medium through which 
the free wind may flow as a supergalactic wind. The energy required to achieve such a 
task, depends strongly on the ISM density distribution. As shown by Silich & Tenorio-Tagle 
(2001) the energy required to burst into the inter-galactic medium out of a fast rotating 
flattened galaxy is orders of magnitude smaller than that required to exceed the dimensions 
of a slow rotating and more spherical density distribution (see Figure 3). 

Given their large UV photon output and mechanical energy input rate, SSCs are now 
believed to be the most powerful negative feedback agents in starburst galaxies, leading not 
only to a large-scale structuring of the ISM and to limit star formation, but also to be the 
agents capable of establishing as in M82 a supergalactic wind, thereby removing processed 
material from galaxies and causing the contamination of the IGM (Tenorio-Tagle et al. 2003). 
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Fig. 3. — Energy estimates. The log of the critical mechanical luminosity, and of the starburst 
mass (right-hand axis), required to eject matter from galaxies with a Mjsm i n the range 10 6 
- 10 9 M Q . The lower limit estimates are shown for galaxies with extreme values of rotation 
(typical of spirals) and for two values of the intergalactic pressure Picu/k = 1 cm -3 K (solid 
lines) and Picu/k = 100 cm" 3 K (dashed lines). The upper limit estimates are for galaxies 
without rotation and thus presenting an almost spherical mass distribution. The resolution 
of our numerical search is AlogL cr = 0.1. Each line should be considered separately as 
they divide the parameter space into two distinct regions: a region of no mass loss (below 
the line) and a region in which blowout and mass ejection occurs (above the line). Also 
indicated on the right-hand axis are the energy input rates required for a remnant to reach 
the outer boundary of a halo with mass 10 8 M and one with mass 10 7 M , for the case of a 
gravitational potential provided by Mjjm — 9.1 xlO 9 M that can hold an Mjsm — 10 9 M . 
The filled squares and triangles represent analytical energy input rate estimates. 
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3.1. The inner structure of M82 

The biconical outflow of M82, the nearest example of a supergalactic wind (SGW), 
displays a collection of kpc long optical filaments embedded into an even more extended 
pool of soft X-ray emission. The outflow is known to extend even further, reaching the 
"H a cap" at 11 kpc from the nucleus of M82 (see Devine & Bally 1999). Both features 
have been partly explained either with the results of Chevalier & Clegg (1985) model of 
an adiabatic, freely expanding, stationary wind and/or by the remnant of a large-scale 
superbubble evolving into the ISM and the halo of the galaxy (see for example Suchkov et 
al. 1994). Both explanations, based on the energetics of a single stellar cluster, fail to explain 
the detailed inner structure of the outflow. The elongated filaments for example, are now 
known to emanate from the central starburst and are not the result of a limb brightened 
superbubble outer structure (Ohyama et al 2002). Note also that the filaments cannot be 
reconcile with instabilities in the large-scale supershell, driven by matter entreinment, which 
in the models occurs at large distances, kpc from the energy source. Furthermore, the 
stationary superwind solution of Chevalier & Clegg leads to a laminar flow and not to gas 
condensation or to a filamentary structure at all. Also, as shown by Strickland & Stevens 
(2000) it has failed to matched the X-ray luminosity of the M82 superwind. Note also that 
the adiabatic assumption, central in the model of Chevalier & Clegg, has recently been 
shown to be inapplicable in the case of massive and concentrated starbursts (Silich et al. 
2003, 2004). Another important issue not accounted by most of the numerical simulations 
is the size of the waist of the biconical structure (150 pc radius in the case of M82), which 
in all calculated cases under the assumption of a single source of energy, (perhaps with the 
only exception of Tenorio-Tagle & Munoz-Tuhon 1997, 1998, which account for the infall of 
matter into the central starburst), also end up with a remnant that presents a wide open 
waist along the galaxy plane (see figures in Tomisaka & Ikeuchi, 1988; Suchkov et al., 1994). 
Further arguments regarding the disagreement between theory and observations are given in 
Strickland & Stevens 2000, Strickland, Ponman & Stevens 1997, and in Tenorio-Tagle et al. 
2003. 



3.2. The physics of supergalactic winds 

So far, all calculations in the literature have assumed that the energy deposition arises 
from a single central cluster that spans several tens of pc, the typical size of a starburst. 
Following however, the indisputable observational findings with HST, we have recently made 
the first attempt to calculate the hydrodynamics that result from the interaction of the winds 
from neighboring young compact clusters present in a galaxy nucleus (see Tenorio-Tagle et 
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al 2003). Several aspects were considered in our two dimensional approach to the problem. 
Among these, the metallicity of the superwind matter was shown to have a profound impact 
on the inner structure of supergalactic winds. Full three dimensional calculations of the 
interaction of multiple SSC winds are now underway. 

Several two dimensional calculations using as initial condition CC85 adiabatic flows have 
been performed with the explicit Eulerian finite difference code described by Tenorio-Tagle 
& Munoz-Tuhon (1997, 1998). This has been adapted to allow for the continuous injection 
of multiple winds (see below). 

We have considered the winds from several identical SSCs, each with a mechanical 
energy deposition rate equal to 10 40 erg s _1 . The energy is dumped at every time step 
within the central 5 pc of each of the sources following the adiabatic solution of Chevalier 
& Clegg (1985). The time dependent calculations do not consider thermal conductivity but 
do account for radiative cooling, with a cooling law (Raymond et al. 1976) scaled to the 
metallicity assumed for every case. 

Figure 4 presents the results for which the assumed metallicity of the winds was set 
equal to 10Z Q , justified by the high metallicity outflows expected from massive bursts of star 
formation (see Figure 1). The winds from the SSCs are exposed to suffer multiple interactions 
with neighboring winds and are also exposed to radiative cooling. For the former, the issue 
is the separation between neighboring sources and for the latter the local values of density, 
temperature and metallicity. Radiative cooling would preferably impact the more powerful 
and more compact sources, leading to cold (T ~ 10 4 K) highly supersonic streams. 

Figure 4 shows three equally powerful (Lsc = 10 40 erg s _1 ) superstellar clusters sitting 
at 0, 60 and 90 pc from the symmetry axis. All of them with an R sc = 5 pc, produce 
almost immediately a stream with a terminal velocity equal to 1000 km s _1 . At t = yr 
the three clusters are embedded in a uniform low density (p = 10~ 26 g cm~ 3 ) medium. Thus 
our calculations do not address the development of a superbubble, nor the phenomenon of 
breakout from a galaxy disk or the halo, into the IGM. The initial condition assumes that 
prior events have evacuated the region surrounding the superstellar clusters, and we have 
centered our attention on the interaction of the supersonic outflows. 

Figure 4 shows the development of a SGW until it reaches dimensions of one kpc, 
together with the final temperature structure splitted into the four temperature regimes: 
The regime of H recombination 10 4 K - 10 5 K, followed by two regimes of soft X-ray emission 
10 5 K - 10 6 K, and 10 6 K - 10 7 K and the hard X-ray emitting gas with temperatures between 
10 7 K - 10 s K. 

The crowding of the isocontours in the figures indicates steep gradient both in density 
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Fig. 4. — Two dimensional superwinds. The first four panels represent cross-sectional cuts 
along the computational grid showing: isodensity contours with a separation A log p = 0.1 
and the velocity field for which the longest arrow represents 10 3 km s _1 . The following four 
panels display isotemperature contours, within the range 10 4 K - 10 5 K, 10 5 K - 10 6 K, 10 6 
K - 10 7 K and 10 7 K - 10 8 K, respectively. Each of the superwinds has a power of 10 41 erg 
s _1 and a radius of 5 pc. The evolution of each wind starts from the adiabatic solution of 
CC85. The plots displays the whole computational grid: 100 pc x 1 kpc. The assumed 
metallicities were Z = 10 Z . The evolutionary times of the first four panels is 1.79 x 10 5 
yr, 4.82 x 10 5 yr, 1.05 x 10 6 yr and 1.39 x 10 6 yr, respectively. 
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or in temperature and velocity, and thus traces the presence of shocks and of rapid cooling 
zones. 

The interaction of neighboring supersonic winds causes the immediate formation of their 
respective reverse shocks, and of a high pressure region right behind them. The pressure 
(and temperature) reaches its largest values at the base of the interaction plane, exactly 
where the reverse shocks are perpendicular to the incoming streams. The high pressure gas 
then streams into lower pressure regions, defining together with radiative cooling, how broad 
or narrow the high pressure zones, behind the reverse shocks, are going to be. 

This also happens if cooling is fast enough, the oblique reverse shocks rapidly acquire 
a standing location, however in these cases, the loss of temperature behind the shocks is 
compensated by gas condensation, leading to narrow, dense and cold filaments. The drastic 
drop in temperature occurs near the base of the outflow, where the gas density is large and 
radiative cooling is exacerbated. The dense structures are then launched at considerable 
speeds (~ several hundreds of km s _1 ) from zones near the plane of the galaxy. These dense 
and cold structures are easy target to the UV radiation produced by the superstellar clusters 
and thus upon cooling and recombination are likely to become photoionized. Note however 
that as the free winds continue to strike upon these structures, even at large distances from 
their origin, the resultant cold filaments give the appearance of being enveloped by soft X-ray 
emitting streams. 

All of these shocks are largely oblique to the incoming streams and thus lead to two major 
effects: a) partial thermalization and b) collimation of the outflow. These effects result from 
the fact that only the component of the original isotropic outflow velocity perpendicular 
to the shocks is thermalized, while the parallel component is fully transmitted and thus 
causes the deflection of the outflow towards the shocks. This leads both, to an efficient 
collimation of the outflow in a general direction perpendicular to the plane of the galaxy, 
and to a substantial soft X-ray emission associated with the dense filamentary structure, 
extending up to large distances (kpc) from the plane of the galaxy. In the figures one can 
clearly appreciate that the oblique shocks, confronting the originally diverging flows, lead to 
distinct regions where the gas acquires very different temperatures, allowing for radiation in 
different energy bands. 

From our results it is clear that a plethora of structure, both in X-rays and in the optical 
line regime, as in M82, may originate from the hydrodynamical interaction of neighboring 
winds. The interaction leads to multiple standing oblique (reverse) shocks and crossing 
shocks able to collimate the outflow away from the plane of the galaxy. In our two dimensional 
simulations, these are surfaces that become oblique to the diverging streams and thus evolve 
into oblique shocks that thermalize only partly the kinetic energy of the winds causing a 
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substantial soft X-ray emission at large distances away from the galaxy plane. Surfaces that 
at the same time act as collimators, redirecting the winds in a direction perpendicular to the 
plane occupied by the collection of SSCs. Radiative cooling behind the oblique shocks leads, 
as soon as it sets in, to condensation of the shocked gas, and thus to the natural development 
of a network of filaments that forms near the base of the outflow, and streams away from the 
plane of the galaxy to reach kpc scales. Under many circumstances the filaments develop 
right at the base of the outflow and for all cases the prediction is that they are highly 
metallic. Hydrodynamic instabilities play also a major role on the filamentary structure. 
Nonlinear thin shell instabilities as studied by Vishniac (1994) as well as Kelvin Helmholtz 
instabilities, broaden, twist and generally shape the filaments as these stream upwards and 
reach kpc scales. 

We thus postulate that if a collection of SSCs is sitting in a preferential plane, most of 
the injected energy would be channeled in a direction perpendicular to the plane of the host 
galaxy. This is achieved naturally as a consequence of the plethora of oblique and crossing 
shocks that redirect the initially isotropic winds. Collimation thus occurs without the need 
of a thick interstellar matter disk or a torus. 

Our considerations point at a new set of possible parameters that profoundly impact 
the development of supergalactic winds. These are: 

• The number and location of superstellar clusters within a galaxy nucleus. 

• The intensity of star formation, or stellar mass in every superstellar cluster, which 
defines their mechanical luminosity, 

• The age of the SSC, which impacts on the metallicity and thus on the cooling of the 
ejected matter. 

All of these are relevant new parameters that may promote, as in M82, the inner struc- 
ture of a supergalactic wind: the co-existense of X-rays and dense filaments, even at large 
distances from the sources of energy. Parameters that may promote self-collimation and with 
it the narrow waist of the biconical outflow. All of these features have been confirmed with 
full 3-D calculations, subject of a forthcoming communication. Note also that our results 
led to the full analysis of the HST data of M82 and to the discovery of 197 stellar clusters 
in the nucleus of M82 (see Melo et al. 2005). 
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4. Positive feedback from massive and compact SSC 

The location of the threshold line in the Esc vs size (Rsc) diagram (see Figure 5), the 
line that defines whether or not a wind is inhibited, depends on several variables. It depends 
on the size of the star-forming region (Rsc) an d the metallicity of the ejected gas, which has 
a strong impact on the cooling curve. It also depends on the assumed E S c/M S c or adiabatic 
terminal speed (i^) of the wind. The latter is also bound to the usual assumption that the 
energy deposited by SN is always 10 51 erg but the mass of the stars exploding within the 
cluster ranges from, say, 100 M Q to 8 M Q and so the injection speed (similar to i>oo) and 
the deposited amount of matter are also functions of time. Another factor that strongly 
affects the location of the threshold line is the thermalization efficiency (e) which defines 
the fraction of the mechanical energy that upon thermalization, can be evenly spread within 
the cluster volume. Estimates of e by several authors lead to values between 1 (Chevalier 
& Clegg; 1985) and 0.03 (see Melioli & Del Pino 2004 and references therein) and depends 
simply on the proximity of the sources undergoing winds and SNe, which through radiation 
may reduce the amount of energy available after thermalization. We have shown that there 
are three different types of solutions: SSCs far away from the threshold line (low mass, low 
energy clusters) undergo a quasi-adiabatic evolution well described by the Chevalier & Clegg 
(1985) solution. More energetic clusters are to have strongly radiative winds. Cooling hardly 
affects their velocity (v w ~ v^) and density distribution (p w oc r~ 2 ), but their temperature 
instead of falling as r -4 / 3 , it falls rapidly to T w ~ 10 4 K close to the SC boundary and the 
more so, the closer they are to the threshold line (see Figure 2). The strongly radiative 
winds around such clusters lead, compared to the adiabatic solution, to very much reduced 
X-ray envelope sizes. The third solution is for clusters above the threshold line. These would 
have their winds inhibited and as shown below, this turns them into very efficient positive 
feedback star-forming agents. 

The facts above the threshold line are that radiative cooling drastically diminishes the 
sound speed esc an d the pressure gradient across the SSC volume, inhibiting the possibility 
of a wind. Radiative cooling upsets then the stationary condition in which the deposited 
matter (Msc) has to equal, at all times, the amount of matter streaming out of the SC 
volume 

M sc = ZEsc/vlc = 4irR sc pscc S c (8) 

where v oo is the resultant wind terminal speed in the absence of radiative cooling. As soon 
as this happens, the mass returned by the stars (Msc) begins to accumulate, promoting 
larger densities and an even faster cooling within the SSC volume. Following this trend, the 
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Fig. 5. — The threshold line. The location of the threshold line, the line that divides the Esc 
vs size (Rsc) diagram into two distinct areas. There the matter deposited within the SSC 
volume, through winds and SNE, ends up either streaming away as a stationary (adiabatic 
or radiative) wind (below the line) or it accumulates to end up steadily being driven into 
new episodes of star formation (above the line). Different locations of the threshold line, for 
an assumed full thermalization efficiency (e = 1), are display for an ejecta metallicity value 
equal to solar and ten times solar. The change in location for different values of e is also 
indicated. 
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accumulated gas density ultimately fulfill the gravitational instability criterion causing the 
collapse into a new stellar generation. 

For clusters above the threshold line, Figure 6 shows how the density of the accumulating 
gas (p ac = 3Msct/4:7iRg C ; where t the evolution time), grows as a function of time within 
the SSC volume, until the moment when the density of the accumulating gas exceeds the 
gravitational instability criteria 

pj ~ 2.3 x 10- 20 (j^j (R scl y 2 g cm-' (9) 

where pj is the Jeans density and Rsci is the SSC radius in pc units. At that moment, when 
Pac = Pj collapse will inevitably proceed. Note that for a given SSC (with a fixed volume), 
P jeans is only a function of temperature and this is set, at least initially, by photoionization. 

For massive clusters the initial ample supply of UV photons exceeds at first the number 
of recombinations within the volume occupied by the reinserted gas and the resultant HII 
region, given the large metallicities of the ejecta, is here assumed to rapidly approach an 
equilibrium temperature T HII ~ a few 10 3 K. At these temperatures the sound speed (< 10 
km s _1 ) remains well below the escape speed and the reinserted gas would inevitably continue 
to accumulate to rapidly (within 1.5 x 10 6 yr) reach the value of the Jeans density for a gas 
at say, Thii = 5000 K, and collapse into a new stellar generation. The event gives rise to a 
new phase of matter accumulation, which once more will rapidly approach pj (for T— 5000 
K) and undergo collapse within a free-fall time, of the order of 10 5 yr, while transforming 
~ 2 x 10 5 M into stars. All stellar generations resultant from mass accumulation within 
the SSC volume have here been assumed to also acquire a Salpeter IMF with similar upper 
and lower mass limits as those imposed to the main superstellar cluster, and their resultant 
properties (mechanical energy and UV photon output) have been added to those produced 
by the main cluster. 

A few, almost identical, stellar generations are expected from the accumulation process 
(solid rising lines in Figure 6), every time that the accumulated gas density p ac reaches 
Pj(5000K) (dashed line in Figure 6). The situation changes slightly when the number 
of ionizing photons (iV°), despite the added contribution of secondary stellar generations, 
becomes insufficient to fully ionize the matter accumulated within the star cluster volume. 
This is due to the evolution of the main cluster, whose UV photon output begins to fall as t~ 5 
after ~ 3.5 Myr. Figure 6 shows pun (thin solid line), the maximum density within the SSC 
volume that can be supported fully ionized by the UV radiation produced by the evolving 
cluster (pun = (3iV° p 2 / (AttR 3 sc (3)) - 5 ; where p = lAmn and (3, the recombination coefficient 
to to all levels above the ground level = 2.59 xl0~ 13 cm~ 3 s _1 ). During the accumulation 
process, once p ac exceeds Phii, the ionized volume begins to shrink to end up as a collection 




Fig. 6. — Positive star-forming feedback. The figure shows the rapid growth of the reinserted 
gas density within the SSC volume (rising solid lines), interrupted as it reaches the gravi- 
tational instability criterion (pj), leading to a new stellar generation and to a new phase of 
matter accumulation. The situation changes slightly as the available photon flux becomes 
unable to sustain all the SSC volume fully ionized, allowing for recombination and further 
cooling of the recombined matter, and thus to a much lower limit of the Jeans instability 
criterion. This allows the accumulation time to become comparable to the free-fall time, 
establishing a new stationary condition in which Msc — SFR. The density value {p gas ) re- 
quired for this condition is also indicated in the figure as well as the maximum density value 
(Phii) that can be supported fully ionized within the SSC volume, with the available UV 
photon flux. 



-18- 



of ultra compact HII regions around the most massive stars left within the cluster, while 
the bulk of the ejected material, now recombined, continues to cool, approaching rapidly 
a temperature ~ 100 K. The characteristic cooling time (£a = 3kT/2An; where A is the 
cooling rate, a function of T and Z) is also very small, ~ 1.5 xlO 5 yr, and is to become 
even shorter as matter continues to accumulate. Matter is at all times uniformly replenished 
within the whole SSC volume, and thus the gas density presents an almost uniform value. 
However, the accumulating gas now has two different temperatures (T HII and 100 K) and 
as p ac grows and the fraction of the ionized volume (fun = (3-/V°/x 2 )/(47ri?| c/ 9p^ c )) shrinks, 
the size of cold condensations (at 100 K) able to become gravitationally unstable and their 
free-fall time also become smaller. 

The drop in the number of ionizing photons and the consequent growth of the neutral 
volume, lead then to a second important condition in which the characteristic accumulation 
time 

Tac = -np gas R 3 sc (l - f HII )/M (10) 
becomes equal to the free-fall time 



This condition defines p gas , the density above the Jeans instability limit for a neutral con- 
densation at 100 K (pj(100K)): 



1/3 / • \ 2/3 

27M 2 \ „ , 18 / M \ (R sc 



and thus once p ac becomes equal to p gas , a new stationary solution becomes possible. 

Everything happens very rapidly, compared to the evolution time-scale of the parental 
cluster (~ 40 Myr), and almost at the same time. The ejected matter is thermalized within 
the SSC volume and immediately begins to cool. At the same time that it accumulates 
making cooling even faster. This allows it to rapidly reach the required p gas value, above the 
Jeans instability limit, that warrants its collapse in a similar time-scale, while the collapsing 
material is replenished by the newly ejected matter. When this happens, the mass deposition 
rate from the cluster becomes equal to the rate of star formation. Gravitational collapse and 
star formation within the star cluster volume and with the matter injected by all sources, 
drive in this way a new stationary condition through a new era of quasi-continuous star 
formation in which M sc is now equal to the star formation rate (SFR). Further details 
describing this phase are given in Tenorio-Tagle et al. 2005. 
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5. Conclusions 

Superstellar clusters are certainly the main mode of massive star formation in starburst 
and interacting galaxies. We have reviewed here how is that they work and the possible 
impact that they may have into the surrounding ISM. We have revised the assumptions of 
Chevalier & Clegg (1985) dealing with thermalization and the flow requirements to estab- 
lish a stationary superwind emanating from these sources. We have also emphasized the 
importance of radiative cooling and how does it affect the superwind X-ray envelopes. 

Calculations in the literature have left clear the fact that single energy sources lead to 
superbubbles and supershells. However, to produce a supergalactic wind with a detailed 
inner structure as in M82, multiple sources, seem to be required. 

We have also shown that a straight forward extrapolation towards the most massive 
and compact coeval clusters is not valid. When radiative cooling becomes significant within 
the SSC volume itself, then instead of driving a superwind able to disperse the surrounding 
ISM and even channel its way into the IGM, events that have make them been regarded as 
negative feedback agents, they become in fact extreme examples of positive star formation 
feedback. 

Massive and compact coeval clusters appear in the E S c vs cluster size diagram above the 
threshold line, in the region where radiative cooling inhibits the development of stationary 
superwinds. In such cases the matter reinserted, through stellar winds and supernovae, is 
unable to escape and after a short phase of matter accumulation, a new stationary solution 
in which M$c becomes equal to the SFR is rapidly met. A positive feedback condition in 
which new stellar generations result in situ, from the collapse of the matter reinserted by 
the sources evolving within the star cluster volume. 

The massive concentrations imply a high efficiency of star formation which permits even 
after long evolutionary times the tight configuration that characterizes them, despite stellar 
evolution and its impact through photo-ionization, winds and supernova explosions, believed 
to efficiently disperse the gas left over from star formation. It is thus the self-gravity that 
results from the high efficiency what keeps the sources bound together. 

The secondary star formation process while causing a faster mass deposition rate, drives 
the SFR to grow from 0.1 to 0.25 M Q yr" 1 over the parent cluster supernova phase (~ 40 
Myr). The continuous reprocessing of the ejected material leads effectively to a continuous 
transformation of the high mass stars into a low mass (< 8 M ) population, keeping constant 
the total mass of the stellar component. 

A central issue, regarding ISM studies, is the fact that the more massive and compact 
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clusters (as those detected by HST and other large ground-based telescopes), are unable to 
generate superwinds and shed their metals into the ISM or the IGM. Their evolution leads 
to many stellar generations and thus to a mixture of stellar populations, all contaminated 
by the products from former stellar generations. An exacerbated episode of star formation 
that leaves no trace of its evolution in the ISM. 
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